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Abstract 

Background: Three gametoclonal plants of Citrus Clementina Hort. ex Tan., cv. Nules, designated ESP, FRA, and ITA 
(derived from three labs in Spain, France, and Italy, respectively), were selected for cytological and molecular 
characterization in order to elucidate genomic rearrangements provoked by haploidization. The study included 
comparisons of their ploidy, homozygosity, genome integrity, and gene dosage, using chromosome counting, flow 
cytometry, SSR marker genotyping, and array-Comparative Genomic Hybridization (array-CGH). 

Results: Chromosome counting and flow cytometry revealed that ESP and FRA were haploid, but ITA was 
tri-haploid. Homozygous patterns, represented by a single peak (allele), were observed among the three plants at 
almost all SSR loci distributed across the entire diploid donor genome. Those few loci with extra peaks visualized as 
output from automated sequencing runs, generally low or ambiguous, might result from amplicons of paralogous 
members at the locus, non-specific sites, or unexpected recombinant alleles. No new alleles were found, suggesting 
the genomes remained stable and intact during gametogenesis and regeneration. The integrity of the haploid 
genome also was supported by array-CGH studies, in which genomic profiles were comparable to the diploid 
control. 

Conclusions: The presence of few gene hybridization abnormalities, corroborated by gene dosage measurements, 
were hypothetical^ due to the segregation of hemizygous alleles and minor genomic rearrangements occurring 
during the haploidization procedure. In conclusion, these plants that are valuable genetic and breeding materials 
contain completely homozygous and essentially intact genomes. 
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Background 

Haploid plants or their derivatives, e.g. doubled haploid 
(DH) or tri-haploid (TH), are valuable in conventional 
breeding and genetic studies. However, most Citrus ge- 
nomes are highly heterozygous, and it is practically impos- 
sible to develop homozygous lines through conventional 
hybridization, due to sexual incompatibility, nucellar 
embryony, severe inbreeding depression, and long juve- 
nility. Gametic embryogenesis is a single-step approach to 
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produce homozygous clones from heterozygous parents 
[1-7], from which most Citrus haploids were generated. 

In situ parthenogenesis induced by irradiated pollen, 
followed by in vitro embryo culture has been reported in 
Citrus [8-11]. Haploid plantlets have been recovered by 
anther culture from Poncirus trifoliata L. Raf. [12] and 
C. madurensis Lour. [13]. A doubled haploid plantlet has 
been obtained from the hybrid No. 14 of C. ichangensis x 
C. reticulata [14]. Haploid plantlets and highly embryo- 
genic haploid calli were recovered from C. Clementina 
Hort. ex Tan. [15-17]. Haploid, but albino embryos, arose 
from cultures of 'Mapo' tangelo (C. deliciosa x C. paradisi) 
[18]. In other reports, haploid and diploid calli, embryos 
and leafy structures, but no green plants, were produced 
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from culture of C. limon L. Burm. f. anthers [19]. Haploid 
embryos of Clausena excavata [20] and homozygous 
short-lived plantlets of Rhode Red Valencia sweet orange 
[21] have also been reported. 

The objective of this work was to elucidate the effect 
of haploidization in the genome structure of three differ- 
ent gametoclonal plants of Citrus Clementina Hort. ex 
Tan., cv. Nules. To compare their genomes, the three 
gametoclones obtained by gynogenesis or by pollen em- 
bryogenesis, were freely provided by research groups in 
Spain (Navarro), France (Ollitrault), and Italy (Germana). 
The tissues and DNA samples from the three candidate 
plants have been analyzed and characterized using vari- 
ous technologies and methods by laboratories in several 
institutions worldwide to assure that they are free of 
large deletions or other defects, as well as to confirm 
their homozygosity (mono-allelic at any locus analyzed). 
Specifically, candidate tree chromosome numbers were 
verified for their ploidy levels. The candidate genomes 
were evaluated using genomic or EST-derived SSR 
markers and microarray technology. The collaborative 
results on the three materials are reported here in 
detail. 



Methods 

Plant material 

Three gametoclonal plants, designated ESP, FRA, and 
ITA respectively acquired in the lab of Navarro (Spain), 
Ollitrault (France), and Germana (Italy), were all derived 
from Citrus Clementina Hort. ex Tan., cv. Nules and pre- 
liminarily shown to be homozygous based on some se- 
lected loci. They were obtained by in situ parthenogenesis 
induced by irradiated pollen followed by in vitro embryo 
culture, or by pollen embryogenesis. Specifically, ESP was 
through in v/vo-induced gynogenesis by pollination of 
Nules Clementine with irradiated pollen of Fortune man- 
darin followed by in vitro embryo rescue [22], FRA also 
through gynogenesis by pollination in the field with irradi- 
ated Meyer lemon (Citrus meyeri Y. Tan.) pollen [11], and 
ITA was obtained through anther culture of C. Clementina 
cv. Nules [15,17]. ESP was previously characterized as a 
haploid [22]. All three plants were much less vigorous 
than the heterozygous mother plant, as revealed by leaf 
size and growth habit (Additional file 1: Figure SI). Sam- 
ples from all three plants were sent to the respective la- 
boratories of the collaborators for the specific analyses to 
which each group had committed. 

Chromosome number determination 

Root- and shoot-tip chromosome counting was conducted 
using DAPI (4,6-diamidino-2-phenylindole) and hematoxylin 
staining techniques [23], respectively. 



SSR genotyping and analysis 

SSR markers presumably or evidently heterozygous in 
the diploid Clementine control were selected for geno- 
typing and analysis to determine if the three plant ge- 
nomes were homozygous and complete at the various 
loci examined. Five laboratories were involved (Table 1), 
and primer sets, amplification conditions, and separation 
methods were summarized. At CIRAD, INRA, and IVIA, 
amplifications were performed according to Froelicher 
et al [24] in a thermocycler (PTC-200, MJ Research) 
using 10 ng of citrus DNA, 0.2 \iM of each primer and 
0.8 unit of Taq polymerase (Goldstar, Eurogentec). The 
annealing temperature for all primers was 55°C. The 39 
EST-SSR primers used in this study were selected based 
on their representation of each of the linkage groups de- 
fined in a Clementine genetic map [25,26]. At the Uni- 
versity of Florida - CREC, 40 EST-SSR markers, likewise 
well distributed among the linkage groups in the diploid 
Clementine genetic map [25], were chosen for genotyp- 
ing. All the genotyping, computing, and scoring proce- 
dures were previously described in detail [27,28]; 
similar methods were used at CCSM [29]. At the Uni- 
versity of California, Riverside (UCR) amplifications 
were performed essentially according to Barkley et al. 
[30], except that PCR products were labeled by adding 
a 19 or 20-base M13 tail to the 5' end of one sequence- 
specific primer and including an M13F or M13R primer 
carrying a dye label (LiCor IRD700 or IRD800) in each 
PCR reaction [31]. 



Array-CGH analysis 

Array-CGH was performed as described in Rios et al 
[32]. Genomic DNA was isolated from leaves [33]. Four 
Cy3- or Cy5-labelled samples from each gametoclo- 
nal plant were co-profiled on four 20 K Citrus cDNA 
microarrays containing 21240 EST, using Cy5- or Cy3- 
labelled control genomic DNA, respectively [34]. To 
prepare labelled probes, Cy3- or Cy5-dCTP fluorescent 



Table 1 Summary of SSR marker analysis of haploids 
from Clementine 



No. of PCR 
products in 
diploid 
Clementine 


No. of PCR 
products in 
ITA, ESP, 
and FRA 


CREC 


INRA 


CCSM 


CIRAD/ 
IVIA 


UCR 


Total 


1 


1 


0 


0 


0 


0 


45 


45 


1 


0 or 1 


0 


0 


0 


0 


2 


2 


2 


1 


40 


39 


9 


41 


57 


186 


2 


1 or 2 


0 


0 


1 


0 


3 


4 


Total markers 




40 


39 


10 


41 


107 


237 


Anomalous 




0 


0 


1 


0 


5 


6 
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nucleotides (Amersham Biosciences) were incorporated 
directly in control and gametoclonal genomic DNA 
(2 (ig) using BioPrime Array CGH Genomic Labelling 
System (Invitrogen). Each pair of purified Cy5 and Cy3 
probes (about 50 \A each) was combined and mixed with 
30 \ig CoM DNA (Invitrogen), 100 \ig yeast tRNA 
(Invitrogen), and 346 \A TE buffer pH 7.4. Samples were 
concentrated with a microcon YM-30 filter (Millipore), 
and SSC buffer and SDS were added to reach a final vol- 
ume of 60 [A containing 3.4x SSC and 0.3% SDS. For 
microarray hybridization, the probe mixture was dena- 
tured by heating at 97°C for 5 minutes, and immediately 
incubated at 37°C during 30 minutes to block repetitive 
DNA sequences. Hybridization mixture was applied to a 
37°C pre-warmed hybrid-slip (Sigma), and a pre-warmed 
array slide was lowered onto the mix. Microarrays were 
hybridized in darkness at 65°C overnight (16-20 hours) 
using a glass array cassette following manufacturers 
instructions (Ambion). To prevent evaporation of hybri- 
dization solution during incubation, 5 \A of 3x SCC were 
poured into the reservoir inside the cassette chamber. Fol- 
lowing hybridization, microarray slides were placed in a 
rack and the cover slip removed by 10 minutes immersion 
in a washing chamber containing 2x SSC and 0.03% SDS 
at room temperature (RT). Microarray slides were passed 
through a series of washes on a shaking platform. Wash 
series were as follow: 2x SSC, 0.03% SDS for 5 min at 
65°C, followed by lx SSC for 5 min at RT, and 3x 15 min 
washes in 0.2 x SSC at RT. Microarray slides were dried by 
centrifugation for 5 min at 300 rpm. Arrays were immedi- 
ately scanned at 5 \xm. Cy3 and Cy5 fluorescence intensity 
was collected by using a ScanArray Gx (Perkin Elmer). 
The resulting images were overlaid and spots identified by 
the ScanArray Express program (Perkin Elmer). Spot qual- 
ity was first measured by the signal-to-background me- 
thod with parameters lower limit (200) and multiplier (2), 
and subsequently confirmed by visual test. The results 
were normalized for labeling and detection efficiencies of 
the two fluorescence dyes, prior to determining differen- 
tial gene expression between haploid and diploid citrus 
samples. Intensities of selected spots were transformed 
into log2 (Cy3/Cy5) and data were normalized by the 
locally weighted linear regression (LOWESS) method. 
Genespring vs 7.3 software (Agilent, http./Avwwagilent. 
com) was used to normalize values for each gene and for 
data analysis. Differentially regulated genes were ranked 
on the basis of signal intensity, normalized ratio, flag value 
and variance across 4 replicate experiments. Filtered genes 
identified to be differentially expressed by haploid/ diploid 
signals lower than 0.3 with a P-value not higher than 0.05 
were considered for subsequent gene dosage measure- 
ments. One-way ANOVA, parametric test without the 
assumption of equal variances was used to define dif- 
ferentially expressed genes. 



Gene dosage measurement 

Quantitative real-time PCR was performed on a LightCycler 
2.0 instrument (Roche), using the LightCycler FastStart 
DNA MasterPLUS SYBR Green I kit (Roche). Reaction 
composition and conditions followed manufacturers in- 
structions. Each individual PCR reaction contained 2 ng of 
genomic DNA from gametoclonal or diploid control [33]. 
Cycling protocol consisted of 10 min at 95°C for pre- 
incubation, then 40 cycles of 10 sec at 95°C for denatur- 
ation, 10 sec at 60°C for annealing and 10-25 sec at 72°C 
for extension. Fluorescent intensity data were acquired 
during the extension time. Specificity of the PCR reaction 
was assessed by the presence of a single peak in the dis- 
sociation curve after the amplification and through size 
estimation of the amplified product. For gene dosage mea- 
surements, the relative quantification-monocolor analysis 
from the LightCycler Software 4.0 package (Roche) was 
used. This program compares the ratio of a target se- 
quence to a reference DNA sequence, i.e. the sequence in 
the gametoclonal sample with the sequence in a diploid 
wild type sample. PCR and normalized calculations were 
repeated in at least three independent samples from each 
genotype, rendering an averaged estimation (± standard 
deviation) of target gene dosage in the haploids. Primers 
for the reference sequence are provided in Table 2. 

Results and discussion 

Chromosome counts 

After DAPI and hematoxylin staining of chromosomes 
in independent labs, 9 were found in ESP and FRA, and 
27 chromosomes in ITA (Figure 1), confirming their 



Table 2 Primer sequences for each gene 


EST name 


Primer sequence 




3 '-CCC A AGCC AG A^TG ATC A AGGGTC-5 ' 


C04035D02 






3 '-TGG ATGTCAC ACC ACTCC AGC AG AT-5 ' 




3'-AGCGCTC™AATCAACCCGTC-5' 


IC0AAA36DF07 






3 '-GG ATACTGCTG ACTG ATG^GC-5 ' 




3 '-TG ATCTCG^TG AGGGTCCTC-5 ' 


IC0AAA34BC06 






3 '-GC A A^CGCC AC^CAGGGTAA-5 ' 




3'-ATGAAGTGTGAGGGTOCGTO-5' 


C34004E09 






3 '-^CGGTC ATGGTC1TC AG AGGT-5 ' 




3 ^IGGTTCA AG AG AGG AGTTGTTG-5 ' 


IC0AAA74CE10 






3'-GATGCAACACATCAGGTGGGAT-5' 




3 '-TGG A^TGC^GGTGC AC ACTG-5 ' 


C06013D07 






3'-GCTGI 1 1 IU ICAACTACAGATCC-5' 




3'-AGTGGGAmGGTGTGGCAAAC-5' 


C08012E04 






3 '-ACCTACTGG A A ATCTG A AG ACC-5 ' 




3 '-COTCCTCATCCAC^C AGG-5 ' 


IC0AAA56DH07 






3 '-CTG AG AC AG A AGCGC A A AC^G-5 ' 



C04035D02 was used as reference gene. 



Germana et al. BMC Plant Biology 201 3, 13:1 29 
http://www.biomedcentral.com/1471-2229/13/129 



Page 4 of 8 




ITA 




Figure 1 Determination of chromosome number by DAPI (A, C) 
and hematoxylin (B, D) staining of chromosomes showing 
counts of 9 in the haploid plants from France (FRA) and 27 in 
the tri-haploid plant from Italy (ITA). Chromosome count for the 
plant from Spain (ESP) was described in Aleza et al. [22]. 



haploidy and tri-haploidy, respectively. They were further 
confirmed by flow cytometry (Additional file 1: Figure S2). 

SSR marker analysis 

A total of 237 SSR markers were selected, in many cases 
from previous mapping exercises, to represent as broad 
and unbiased coverage of the citrus genome as possible, 
and plant materials were genotyped (Table 1). No SSR 
alleles were detected in the gametoclones that were not 
present in diploid Clementine. At 232 loci the three 
gametoclones had one SSR allele also found in diploid 
Clementine. The gametoclones had the same allele as 
diploid Clementine at 45 of the 47 loci tested at which 
Clementine appeared homozygous. At two "anomalous" 
loci, the Clementine allele was observed in one or more 
of the gametoclones, and no allele was observed in the 
others. These two loci segregated for a null (no amplifi- 
cation) allele in a Clementine hybrid population, so these 
markers are also consistent with all gametoclones having 
complete, homozygous genomes. The gametoclones con- 
tained one of the two Clementine alleles at 183 of the 
187 loci that were heterozygous in diploid Clementine. 
For three loci, the same two PCR products amplified 
from diploid Clementine were also observed in one or 
more of the gametoclones. Segregation of one of these 
loci was studied in a Clementine hybrid population and 



it was shown that these two PCR products segregated as 
a single Mendelian unit. This pattern could be caused by 
a tandem duplication of the amplified region, or by 
annealing of one PCR primer to nearly adjacent sites in 
the template DNA. Segregation of the other two loci has 
not been examined, but they could be explained possibly 
by a similar mechanism. Only one SSR locus revealed 
anomalous results in FRA, while the remainder revealed 
only a single allele product at all other loci surveyed. 

Array-CGH experiment 

In order to detect putative genomic deficiencies, gen- 
omic DNA from ITA, ESP and FRA genotypes and the 
diploid wild type was labelled and hybridized to a 20 K 
citrus cDNA microarray [34] by array-CGH, a procedure 
that previously proved to be useful for the structural 
prediction of large genomic deletions in Citrus 
Clementina [32]. Those ESTs showing a haploid/diploid 
signal ratio of 0.3 or lower with a maximum P-value of 
0.05 for any of the genotypes were selected for further 
analysis. From 13 ESTs fulfilling these conditions, two of 
them were annotated as putative Cu/Zn-superoxide 
dismutase copper chaperones, which are the 5 'end 
(C34004E09) and 3 'end (KN0AAA2CB01) of the same 
citrus cDNA (Table 3). Under-represented ESTs were 
found in the three gametoclonal genotypes in a number 
ranging from 4 in ESP to 8 in ITA with 4 of them jointly 
found in two different genotypes. 

Gene dosage experiment 

Seven of these ESTs were chosen for gene dosage evalu- 
ation by real-time PCR. Gene dosage measurements 
performed in this way confirmed the hybridization data 
presented in Table 3, except in two cases (Figure 2). 
Real-time PCR quantification failed to amplify ESTs 
C06013D07 and IC0AAA56DH07 that exhibited, re- 
spectively, microarray hybridization signal ratios of 0.23 
and 0.30 in ITA, indicating that microarray data of these 
two ESTs was most likely affected by nonspecific cross- 
hybridization. 

Since a gene dosage around one is expected for those 
genes that are neither enriched nor reduced upon 
haploidization, the estimation of gene dosages lower and 
higher than one in certain ESTs argues for the occurrence 
of genomic rearrangements during the haploidization pro- 
cedure, or alternatively the segregation of hemizygotic 
genes, that are present exclusively in one of the two alleles. 
A schematic overview of these and other genomic expla- 
nations for the observed deviations in gene dosage is 
presented in Figure 3. In this Figure, the standard case of 
genes not affected by the haploidization procedure, with 
both haploids having identical alleles and therefore a con- 
stant gene dosage around one is exemplified in panel A. 
In the hemizygotic model, a gene or DNA fragment is 



Table 3 Underrepresented ESTs after array-CGH of haploid genomic DNA 



EST 


Accession 


Genotype ITA 




Genotype FRA 




Genotype ESP 




Description 


Normalized ratio (*) P-value (**) 


Normalized ratio (*) P-value (**) 


Normalized ratio (*) P-value (**) 


IC0AAA36DF07 


DY274520 


2.08(0.98-3.45) 




0.04(0.03-0.06) 


1 .8E-04 


0.04(0.01-0.09) 


5.0E-02 


No annotation available 


IC0AAA34BC06 


DY273568 


0.34(0.26-0.63) 


1 .4E-02 


0.24(0.17-0.31) 


1.5E-03 


0.31(0.28-0.40) 


1 .OE-02 


No annotation available 


C34004E09 


FC930172 


0.95(0.87-1.05) 




0.13(0.09-0.20) 


1 .7E-03 


0.98(0.92-1.05) 




Cu/Zn-superoxide dismutase copper chaperone 


KN0AAA2CB01 


DY256847 


1.00(0.64-1.92) 




0.07(0.01-0.35) 


3.4E-02 


0.78(0.39-2.19) 




Cu/Zn-superoxide dismutase copper chaperone 


C31204A07 


FC874254 


1.04(0.85-1.20) 




0.21(0.10-0.40) 


1 .2E-02 


2.33(1.78-3.87) 




Thaumatin-like protein isoform 3 


IC0AAA74CE10 


DY290159 


0.22(0.18-0.30) 


1.3E-03 


1.65(1.30-2.07) 


1 .4E-02 


0.21(0.18-0.27) 


6.3E-03 


Nematode resistance-like protein 


KN0AAK2BA02 


DY259689 


0.14(0.07-0.19) 


3.6E-03 


1.74(1.24-2.59) 


3.9E-02 


0.21(0.14-0.26) 


1 .4E-02 


No annotation available 


C06013D07 


CX298593 


0.23(0.11-0.40) 


1.1E-02 


1.26(1.11-1.42) 


3.5E-02 


1.37(1.29-1.45) 


1.1E-02 


No annotation available 


C31810G10 


FC924855 


0.18(0.10-0.50) 


1 .7E-02 


1.81(1.35-2.64) 


3.8E-02 


1.35(0.89-2.68) 




Copine l-like protein 


C08012E04 


CX301690 


0.28(0.17-0.54) 


1 .8E-02 


1.87(1.40-2.58) 


2.3E-02 


1.28(0.83-2.99) 




No annotation available 


IC0AAA56DH07 


DY282523 


0.30(0.12-0.63) 


4.0E-02 


1.66(1.43-2.13) 


1 .2E-02 


1.72(1.17-2.13) 




No annotation available 


IC0AAA89CG08 


DY295961 


0.27(0.14-0.45) 


1 .6E-02 


0.43(0.28-0.99) 




1.52(0.87-2.42) 




No annotation available 


KN0AAB2CE09 


DY258014 


0.07(0.03-0.14) 


7.0E-03 


1.59(1.16-2.14) 


4.0E-02 


0.21(0.11-0.66) 




SJCHGC09076 protein 



(*) Normalized haploid to diploid signal ratio. In bold, signal ratios lower than 0.3. Replicates with lowest and highest values are shown in brackets. 
(**) Only Pvalues lower or equal to 0.05 are shown. 
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Figure 2 Relative gene dosage of selected ESTs in the three 
haploids by real-time PCR. Values are average and standard 
deviation of three determinations. Gene dosage is relative to the 
value of the diploid variety. 



absent in one of the alleles leading to gene dosages equal 
to zero in haploids carrying the deleted allele. In the dip- 
loid parental, the PCR amplification of hemizygotic genes 
is originated from only one allele, as the second allele is 
absent (Figure 3B). The relative enrichment of the gene in 
haploid individuals inheriting the full allele causes a two- 
fold increase in gene dosage, whereas haploids carrying 
the null allele have a value close to zero. In this work, 
IC0AAA36DF07, C06013D07 and IC0AAA56DH07 showed 
such a hemizygous-like behavior. In the deletion model, 
as outlined in Figure 3C, the gene is deleted during 
haploidization procedures producing a null allele that is 
not present in the diploid. EST C34004E09, for example, 
was not detectable in FRA but its content in ITA and 
ESP individuals was close to one. These data could be 
explained by a genomic deletion mechanism occurring 
during the haploidization process. Under this model, 
haploids losing the gene during the haploidization do not 
show PCR amplification signal, but those haploids inhe- 
riting an intact allele show a relative gene dosage similar to 
the original diploid. Finally, the remaining three ESTs ana- 
lyzed by real-time PCR show low gene dosage values higher 
than zero in the three haploids (IC0AAA34BC06) or in at 
least one of them (IC0AAA74CE10 and C08012E04). 
Three different structural models were postulated to ex- 
plain these observations. In the polymorphic tandem re- 
peats model (Figure 3D), a tandemly-repeated gene found 
respectively Y and y times in the two alleles, show a gene 
dosage value responding to the equations (2-x)/(x + y) and 
(2-y)/(x + y) in the two alternative haploids. Thus, a 5-fold 
ratio in the number of repeats would produce relative 
gene dosages of 10/6 and 2/6 in the resultant haploids, 
certainly similar to the observed values. Alternatively, 
polymorphic variations in the primer binding sequence on 
the gene might cause allele-specific modifications of PCR 











A -> «- 

— • 1 




0 




D > < . < 2x 


x+y 






Figure 3 Proposed models for the genomic mechanisms 
resulting in lower gene representation in the gametoclones. 

Genes are depicted as blue and green cylinders and primers as red 
(specific) and orange (non-specific) arrows. Yellow flashes label 
deletion limits and the white arrow indicates conversion of standard 
to null allele. In each panel up and down drawings represent two 
alleles of a given locus in the diploid variety. The number or 
equation on the right is the relative gene dosage estimated for a 
haploid line inheriting the corresponding allele located on the same 
level. Gene dosages are always one for the diploid line. For further 
explanations see the text. (A) Standard case; (B) Hemizygotic model; 
(C) Technical deletion model; (D) Polymorphic tandem repeats 
model; (E) Polymorphic sequence model; (F) Hemizygotic model 
and cross-reaction. 



efficiency leading to variable gene dosages (Figure 3E, 
polymorphic sequence model). Another source of mis- 
estimation of gene dosage might result from the combin- 
ation of hemizygosis and non-specific cross-reaction of 
the primers, as presented in Figure 3F that originate al- 
tered determinations of gene dosage in the resulting 
haploid genotypes. Thus, the results confirm that ITA, 
ESP and FRA genomes do not carry important fragment 
deletions or rearrangements, and the few genomic dif- 
ferences observed between the diploid and haploid ge- 
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notypes can be explained to a large extent by the natural 
heterozygosity of the diploid parental 

Conclusion 

In this study, chromosome counting confirmed that ESP 
and FRA were haploid and ITA tri-haploid. Among a total 
of 237 SSR markers, most were selected from previous 
mapping exercises and represented broad and unbiased 
coverage of the citrus genome. 231 markers detected a 
single allele in ITA, ESP and FRA; each allele also existed 
in the diploid Clementine genome. Of the six SSR loci 
with anomalous results, segregation in Clementine was 
studied for three loci and in these cases the anomalous re- 
sults in the haploids were shown to be caused by similar 
anomalies in Clementine. The array-CGH experiment re- 
vealed that only 13 cDNAs had anomalous results among 
more than 20,000 cDNAs on the array. After real-time 
PCR of 7 of these genes, only four showed a gene dosage 
close to zero in one or two candidates, so no relevant gene 
loss was detected in any of the three genomes. Conse- 
quently array-CGH, in addition to all other characteri- 
zation methods employed, provided compelling evidence 
that haploidization of citrus through in situ partheno- 
genesis induced by irradiated pollen followed by in vitro 
embryo culture, or by pollen embryogenesis, does not 
generate substantial genome rearrangements. Therefore, 
these three gametoclones can be used, with no concerns 
regarding their genomic integrity, for genetic studies as 
well as for citrus improvement, for example, through di- 
haploidization. In addition, it is noteworthy that the con- 
clusions reached in this study, that haploidization does 
not disrupt the natural citrus genome structure, provided 
the major basis for the selection of the target citrus gen- 
ome for producing the reference sequence for the inter- 
national citrus research community [22]. 
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cytometry analyses of DNA content for the plants from France (FRA) and 
Italy (ITA). 
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